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Abstract. The results of a recent next-to-leading order QCD analysis of the world data on polarized 
deep inelastic scattering are reported. New parameterizations are derived for the quark and gluon 
distributions, accounting for the massive Wilson coefficient for the charm quarks, and the value of 

as(Mj ) is determined with correlated errors. We obtain a]^'"*^(M|) = 0. 1132 ]'^q Q095 ■ Limits on 

potential higher twist contributions to the structure function gi (x, Q^) are derived. We also compare 
to the results obtained by other groups. 
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INTRODUCTION 

The composition of the nucleoli's spin and the short distance behaviour of the partons 
inside strongly polarized nucleons in general constitutes one of the central research top- 
ics in QCD, with emphasis on the twist-2 contributions. They are explored both with 
perturbative and non-perturbative methods. During the last years the polarized deep- 
inelastic scattering data have further improved [1-15]. In this note we report on a new 
next-to-leading order (NLO) QCD analysis of these data [16], updating earlier investi- 
gations [17]. At large enough four-momentum transfer = —q^, the structure function 
g\{x, Q^) mainly receives twist-2 contributions^ and is related to the polarized twist-2 
parton distribution functions (PDFs). We analyze the structure function ^1 (jc, Q^), which 
is derived from the longitudinal polarization asymmetry accounting for a data-based de- 
scription of the denominator function [19] and corresponding parameterizations for the 
longitudinal structure function, cf. [16]. In the present analysis we include the 0{as) 
contribution due to charm quarks, [20, 21] ^. The structure function g2{x,Q^) is de- 
scribed at leading twist by the Wandzura-Wilczek relation [23, 21]. The parameters of 
the polarized parton densities, which can be measured using the above data sets, are 

mandatorily determined with correlated errors along with the QCD scale ^qq^) ■ 
also analyze potential contributions of higher twist and derive corresponding limits. Fi- 
nally, a phenomenological parameterization of the polarized NLO PDFs is provided in 
terms of grids for the central values and the correlated errors, [24]. 



1 DESY 10-244, SFB/CPP-10-133. Supported in part by DFG SFB-TR-9. 
^ Twist-3 contributions are connected by target mass effects, cf. [18]. 

^ The (9(a^) heavy flavor corrections are only known in the asymptotic region Q^->n?, [22]. 



THE ANALYSIS 



The NLO QCD analysis of the structure function gi (x, Q^) is performed in Mellin space 
following the standard formalism, cf. e.g. [25], including the heavy quark corrections 
[26]. In this representation the evolution equations can be solved analytically, in both 
a fast and numerically precise way. Only one numerical integral around the singularities 
of the solution in the complex plane, located at the real axis left to an upper bound, 
has to be performed to represent gi (x, Q^). As the data are located at low values of 
target mass corrections are applied, cf. [18, 28]. For the deuteron targets a wave function 
correction is performed [29]. The parton distributions at the starting scale Qq=4 GeV^ 
are parameterized by 



with r]i the first moments. The present analysis parameterizes the sea quarks assuming 
approximate flavor SU{3) symmetry. The deep-inelastic data alone cannot resolve the 
flavor dependence of the sea. Taking into account semi-inclusive data [30], and later 
on polarized Drell-Yan and di-muon data, will allow the determination of polarized sea 
quark distributions, similar to the unpolarized case. ^ 77„^ and t]^^ are fixed due to the 
neutron and hyperon-j8 decay parameters F and D, which are very well measured : 



The parameters in (1) cannot all be measured using the present data since for some the 
;^^-fit yields errors larger than 100%. In case of the sea-quark and gluon density Yi is 
found to be compatible with zero. Furthermore the a,-parameters of the distributions A^^ 
and AG are related by about a\G = ^A^^ + 1> which we use. The parameters Juy and Juy 
are fitted in an intial run and are then kept fixed as model parameters. For the large-x 
parameters bAq, and Z?ag we used the relation ^a?s/^ag(Po1) = ^A?s/^AG(unpol) = 1.44 
and determine b^c = 5.61 and b^q, = 8.08 in the fit. In the final fit 8 parameters are 

Nr=4 

determined including Aq'j-.j^ . In Figure 1 we show the four distributions Auy^Ady^Aqg 
and AG at the input scale and compare them to other determinations. In Ref. [16] we 
provide the correlated errors. Due to this one may perform Gaussian error propagation 
for all observales based on polarized parton denstities predicting the PDF errors of these 
quantites. There we also compute a series of moments for the different parton densities, 
see Ref. [16], which can be compared to upcoming lattice simulations. 
The nucleon spin is given by the relation 



xA/;(x,(2o) = T]/A;x«Kl-x)*'(l + T^x) , 



(1) 



riu,.-Vd,= F + D and nu,. + 'nd^,= 3F + D, 
0.928±0.014 and 77^^= -0.342±0.018 . 



(2) 
(3) 



1 



^{AL{x))o + {AG{x))o + Lq+L, 



(4) 



2 



We refrain from carrying out small-x resummations, since yet unknown subleading terms are very likely 
to cancel the leading order effects, cf. [27]. 

^ For first analyses accounting for the flavor dependence of the sea quarks see [31]. 
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FIGURE 1. NLO polarized parton distributions at the input scale = 4.0 GeV^ (solid line) 
compared to results obtained by GRSV (dashed-dotted line) [32], AAC (dashed line) [33], and 
LSS (long dashed line) [34]. The shaded areas represent the fully correlated la error bands 
calculated by Gaussian error propagation. The dotted line indicates the positivity bound using 
the parameterization [35]; from Ref. [16]. 

to the first moments of the polarized flavor singlet and gluon distributions and the quark 
and gluon angular momenta Lq^g. In the present analysis we obtain 

(AE(x))o = 0.216 ±0.079 (5) 
{AG{x))o = 0.462 ±0.430, (6) 

which saturates the required value even for vanishing values for Lq and Lg. However, 
the error on the gluon density is still rather large. Using the grids [24] one may perform 
predictions for polarized observables at hadron colliders which depend on the twist-2 
parton distributions for RHIC or other machines planned for the future, as e.g. for the 
Drell-Yan process or pseudoscalar Higgs boson production [36] ^. 

In deep-inelastic QCD analyses it is important to determine the PDF-parameters at the 
initial scale Qq together with the QCD scale Aqcd since there are strong correlations, 
e.g. between the gluon-normalization and ttsiM^), but also to other parameters. We ob- 



For recent predictions in case of unpolarized /?(/?)/?-scattering at NNLO see [37]. 
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FIGURE 2. A summary of the current measurements of (M|) from unpolarized and polarized 
DIS data, of. Ref. [16] for details. Due to the size of errors we include only the results of 
NNLO and N^LO analyses in the unpolarized case, while those in the polarized case stem from 
NLO analyses. The yellow band marks the weighted world average of 2009, to which, however, 
measurements of as{M^) at NLO, NNLO and N^LO contribute [38]. 
tain 



A, 



(4) 
QCD 



243 ±62 (exp) 



-37 
+21 



(FS) 



+46 
-87 



(RS) MeV 



(7) 



The renormalization (RS) and factorization scales (FS) were varied by a factor of 2 
around Q^. Here we excluded values nj^<l GeV^, unlike in Ref. [17], since at scales 

lower than 1 GeV^ the perturbative description cannot be considered reliable anymore. 
Correspondingly, one obtains 



a, (m|) =0.1132 



+0.0043 
-0.0051 



(exp) 



-0.0029 
+0.0015 



(FS) 



+0.0032 
-0.0075 



(RS). 



(8) 



The errors are much larger than in the unpolarized case, at NNLO, where an accuracy of 
0( 1 %) is reached. Still the central value is lower than the current world average and well 
comparable to the unpolarized values. In Figure 2 we summarize the current status of 
05.5 (M|) measurements in deep-inelastic scattering. We would like to mention the results 
of the unpolarized NS-analysis [39] at N-^LO : 



.7 



a,(M| 



0.1141 



+0.0020 
-0.0022 



and recent combined NS and singlet NNLO analyses [40, 41] 

a^(M|) = 0.1124 + 0.0020 
a,(M|) = 0.1135 + 0.0014. 



(9) 



(10) 
(11) 



The corresponding NLO and NNLO values are a,{M^) = 0.1 148 ± 0.0019 and 0.1 134 ± 0.0020, re- 
spectively, pointing to the fact that the NNLO values are systematically smaller than the NLO values. 
Therefore it is problematic to average a., values extracted at different orders. 
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FIGURE 3. The additive higher twist coefficients Cp {x),Cd{x) and C„ (x) as a function of x. 



Very recently, due to the inclusion of the combined Hl+ZEUS data the latter value 
receives a slight change to 

0.1147 ±0.0012, (12) 



a,(M|) 



reaching now the accuracy of 1 %, [42]. 

We also fitted additive higher twist terms to gi{x,Q^), allowing for a model term 
+C{x) /{Q^ /GeV^) besides the twist-2 contributions, to explore the corresponding 
structures in the region jc < 0.6 for the proton- and deuteron targets, cf. Figure 3. While 
in case of the deuteron target the result is fully compatible with zero, an effect of up 
to 2 a is observed for one out of five bins in case for the proton target. This result is 
indicative mainly, since a measurement of the higher twist contributions would require 
a clear separation of the leading twist terms. This is possible in the unpolarized case, as 
has been shown in Refs. [43, 44]. This also requires even higher order corrections for 
the leading twist contributions. A comparable analysis in the polarized case has to be 
based on much more precise data in a far wider range of which can be obtained at 
future colliders such as the EIC. 
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